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Perfect Lattice Actions for the Gross-Neveu Model 

Wolfgang Bietenholz^, Erich Focht'^ and Uwe-Jens Wiese^ * 

'^Center for Theoretical Physics, Massachusetts Institute of Technology (MIT) 
Cambridge, Massachusetts 02139, USA 

''Institute for Theoretical Physics E, RWTH- Aachen, 52074 Aachen, Germany 
HLRZ c/o KFA JiiHch, 52425 Jiihch, Germany 

We apply the method of Hasenfratz and Niedermayer to analytically construct perfect lattice actions for the 
Gross-Neveu model. In the large A^ limit these actions display an exactly perfect scaling, i.e. cut-off artifacts 
are completely eliminated even at arbitrarily short correlation length. Also the energy spectrum coincides with 
the spectrum in the continuum and continuous translation and rotation symmetries are restored in physical 
observables. This is the first (analytic) construction of an exactly perfect lattice action at finite correlation 
length. 
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1. Introduction 

The worst systematic errors in lattice simula- 
tions are due to cut-off artifacts. For fermionic 
systems they are typically of 0(a), where a is the 
lattice spacing. For a long time, the attempts 
to improve the lattice action such that the arti- 
facts are dampened, yielded only a limited suc- 
cess. Last year, P. Hasenfratz and F. Nieder- 
mayer proposed a new method for approximating 
perfect actions (actions free of cut-off effects) for 
asymptotically free theories [y . It might have the 
potential for a break-through which enables us to 
solve lattice QCD. 

First we note that the fixed point actions 
(FPAs) - obtained from iterated block spin renor- 
malization group transformations (RGTs) on a 
critical surface (see e.g. Q) - are perfect actions. 
However, since the correlation length ^ diverges 
there, they cannot be used in simulations. Kogut 
and Wilson |^ noticed that perfect actions exist 
also off the critical surface. They lie on 'renormal- 
ized trajectories' (RTs) in parameter space, which 
cross the critical surface in a fixed point. The di- 
rection of the tangent there corresponds to a rel- 
evant (or marginal) direction. If we move slightly 



*Based on two talks presented by E.F. and W.B. Work 
supported by Deutsche Forschungsgemeinschaft (D.F.G.) 
and U.S. Department of Energy (D.O.E.) under coopera- 
tive research agreement DE-FC02-94ER40818. 



away from a fixed point in a relevant direction and 
perform RGTs, then we follow the RT. Thus all 
points on it are related to the vicinity of the fixed 
point without occurrence of irrelevant terms. 

For asymptotically free theories, there is only 
one (weakly) relevant (to leading order marginal) 
direction. Hasenfratz and Niedermayer propose 
to follow this direction down to a ^ short enough 
for simulation, i.e. to approximate the RT lin- 
early. The points on this tangent are referred 
to as 'perfect classical actions' and there is hope 
that they still have very small cut-off artifacts. 
The virtue of asymptotically free theories (in this 
context) is that their FPAs can be determined 
classically. 

The pilot project of Hasenfratz and Nieder- 
mayer dealt with the 0(3) nonlinear a model in 
d = 2 and was a full success. Even on a lattice as 
small as 5 X 5 cut-off effects were numerically not 
visible any more. 

The final goal is of course QCD. A collabora- 
tion is presently studying the pure SU{3) gauge 
theory in this respect Q. 

Here we present results for the 2d Gross-Neveu 
model Q , which turned out to be a good testing 
ground for the method when applied to fermionic 
systems. 



2. Fixed point actions for free lattice 
fermions 

First, one of us found a line of FPAs for free 
Wilson fermions H which also applies to a num- 
ber of other lattice fermions (SLAC, Rebbi, etc., 
see 0). However, this set does not include the 
staggered fermions which we are going to use here. 
Their main advantage is the remnant U{1)(E)U{1) 
symmetry which persists at free, local fixed points 
g .g We perform RGTs according to the prescrip- 
tion in yi i.e. we attach pseudo flavors to the cor- 
ners of (disjoint) 2x2 blocks and build block spins 
such that the pseudo flavors don't mix. Therefore 
an odd blocking factor is needed; we chose it to be 

3. Then we have intersecting 9x9 blocks, whose 
centers form a lattice of spacing 3 and which are 
occupied by alternating pseudo flavors. Every 
point in a block with the same pseudo flavor as 
the center contributes to its block spin which lives 
again on the block center. Then we rescale coor- 
dinates and fields. 

We take a lattice of spacing 1/2. For technical 
reasons we initially put the ferniionic pseudo fla- 
vors x^, . . . X*, X^ • ■ • X^ (one-component Grass- 
mann variables) on the centers x of disjoint 2x2 
blocks (spacing 1). Later the fermionic variables 
are moved to the corners (in H we proceeded dif- 
ferently). For the action after a number of RGTs 
we make the general ansatz 



(1) 



S is invariant against translations by integers (but 
not half integers) along the coordinate axes. The 
matrix for the fermions shifted to the corners is 
(in momentum space) denoted by p{k). Its form 
is restricted by charge conjugation, U{1) O U{1) 
symmetry and half integer translations. Its struc- 
ture, resp. the one of its inverse a{k) = p~^(k) 
is 



•(2) 



^For Wilson fermions the full chiral symmetry may be pre- 
served, but then the FPA becomes non local, in agreement 
with the Nielsen-Ninomiya theorem. 
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Starting from the standard action with a^{k) = 
ik^/k'^ (in lattice notation: k^ = 2sin(fc^/2)), we 
arrive at the fixed point given by 
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Here a renormalization parameter in the block- 
ing step had to be chosen in agreement with the 
dimension of x- The case c ^ oo corresponds 
to a (5-function blocking RGT, but we may add 
a 'smearing term' suppressed by 1/c, which does 
not break any symmetry and which helps to op- 
timize the locality, see ||]. 

3. The Gross— Neveu model 

Our main interest was to gain experience with 
interacting fermionic models, the simplest asymp- 
totically free one being the Gross-Neveu model in 
two dimensions (GN). It is a model with 4- Fermi 
interaction which can be linearized by introducing 
a real valued auxiliary field $. With G being the 
usual 4-Fermi coupling the continuum Euclidean 
action for N flavors is 

The 0{2N) x Z{2) symmetry of the model is bro- 
ken spontaneously to 0{2N) and a fermion mass 
mj is generated dynamically. 

In a sensible lattice formulation of the GN 
model with staggered fermions and Yukawa cou- 
pling, the scalars live on the centers of the 
fermionic lattice, i.e. they couple equally to the 
4 pseudo flavors at one plaquette and their spac- 
ing is 1/2. To keep a Gaussian form we do not 
integrate out <1> but block it, simultaneous to the 
fermionic blocking. For this we use the same pat- 
tern: we block them as if they also had pseudo 
flavors, which must not be mixed, see fig.l. 

4. Small field approximation 

For the moment we only include the first order 
in $ and let A^ = 2. Then we parametrize a 
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Figure 1. Geometry of the scalar (crosses) 
and fermionic variables (different symbols for the 
pseudo flavors). The nine auxiliary field variables 
marked with a small circle contribute to the block 
variable in the center of the box. 

general ansatz for the action as: 

^^[x,x,*] = ^{EE^^/'''(^-y)4 (4) 



G 



G' 



+ 2 ^ *'+ ^ ^xWi^ -z,y~ z)xl'^,}. 

where z runs over the cell centers (spacing 1/2). 
For the standard action in momentum space, 
shifting the fermions to the block corners takes 
a to the form 
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The fixed point condition forces us to apply a 5 
blocking in the scalar sector. Doing so we arrive 
at a fixed point, where a still has the same struc- 
ture and its elements become: 
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The 47r antiperiodicity of these elements agrees 
with the central positions of the scalars. Thus 
we have identified the direction of the RT out 
of the critical surface, but only in the parameter 
subspace which is accessible in the framework of 
the small field approximation. 

5. The large N limit 

In general one can hardly go beyond the small 
field approximation to determine the functional 
S*\xi X? '^] analytically. An exception is the limit 
N ^ 00. There the scalar field freezes to a con- 
stant field $0 and we can calculate S* [x, X, '^o] 
nonperturbatively. We parametrize a general 
ansatz as 



5'[x,X,«'o] = 
1 
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where B is the Brillouin zone and V is the space- 
time volume factor (actually infinite). For the 
standard action p is the same as for free fermions 
and A(fc) = $0- Wedefine(5(fc)-|-/3(fc)l = [p{k) + 
A(fc)l]^^ (again the tilde corresponds to fermions 
shifted to the plaquette corners; for A and /3 there 
is no change). Still p and a keep the structure 
given in (H) and we obtain at the fixed point 
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Here the small field approximation allows to in- 
clude small non-zero modes of $(fc) and may be 
used if one is interested e.g. in 1/N effects. The 
result to the first order in $(fc 7^ 0) is given in 



H . The occurring matrix a* {p, q) is more gen- 
eral than the one given in (0) , because we expand 
around any $o now, whereas (ph describes a par- 
ticular expansion around $o = 0. 

6. The perfect operator for the chiral con- 
densate 

Let's perturb S* a little with some operator X: 

Sj[x,X,^ = S*[x,X,^]+jX[x,X:'^] , (9) 

where j <^ 1. We call X a perfect operator and 
denote it by X* if it is an eigenfunctional of the 
RGT to 0{j) 



s'[x\x,'^'] = s*[x',x\'^'] 



(10) 



+nx*[x',x',^'] + 0{f). 

The operator is relevant if 7 > 1. This is the case 
for the chiral condensate (with standard operator 
"llxXxXx)- for dimensional reasons 7 coincides 
with the block factor. 

We specialize again to A^ -^ cx), $ ^ $0 and 
also parametrize X as; 

X[X, X, ^0] =7K-v? fdkxi-kMk) + iyik)l]x{k), 



where jj, has the structure of p and a. From 
the above condition (|l^) for X* we obtain (again 
shifting the fermions to the block corners) 
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(Omitting O(<l>o) we observe that X* is very lo- 
cal.) 



The condition for V^ff to be minimal yields the 
gap equation 
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A4 = —T^[p + Al] being the fermion matrix for 
one set of staggered fermions. For the standard 
action the gap equation takes the form 

9 (2^)2 Jb fc2 
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obtained and discussed before in [y. For the 
FPA it reads 
1 
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First we investigate asymptotic scaling of $0 , see 
fig. 2. Asymptotic scaling corresponds to linear 
curves. We see that asymptotic scaling for the 
FPA is clearly better, but still not what one would 
call 'perfect'. As an alternative consideration of 
asymptotic scaling we compare in fig. 3 the f3 func- 
tions referring to the fermion mass, which we will 
derive below. It shows the same feature. (Asymp- 
totic scaling is given by: amf oc e~'^'^ ^ f3 = 
-jTifdrufg = -.g/TT^.) At ^ = 1 the FPA's (3 func- 
tion already agrees with this up to 5 percent, but 
its behavior is still not perfect. 



standard action 
perfect action 



7. The gap equation 

Now we forget about all small field approxima- 
tions, let TV —> 00 keeping g — GN fixed and ex- 
ploit our exact result (g) . We make use of the di- 
mensional transmutation due to the spontaneous 
symmetry breaking U{N/2) Z(2) -^ U{N/2), 
which gives rise to a fermion mass rrif . We intro- 
duce an effective potential V^ff by 
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Figure 2. Asymptotic scaling behavior of^Q{g). 
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Figure 3. P{g)l g^ versus 1/g. 

However, this doesn't mean very much. 
Asymptotic scaling is not really physical since it 
involves the bare coupling (at cut-off scale). Ac- 
tually meaningful are dimensionless ratios. 

8. The scaling behavior 

Such a ratio shall be constructed now to discuss 
the scaling behavior of the standard action, the 
small field approximation, and finally of the FPA. 
First we derive mf from the exponential decay of 
the correlation function. We put the system in 
a spatial volume L and consider X2 as Euclidean 
time. The correlation function becomes 
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For the standard action this turns into 
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At the pole of the integrand we find the energy 
spectrum: 
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We observe cut-off artifacts of 0{a^). Wc identify 
■nif = E{Q) such that 
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Now we consider the chiral condensate 
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For the standard action the gap equation (|lj) 
leads to (xx) = — $0. Therefore 
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In the continuum limit, this ratio must become 
constant, and indeed it does, since: ^ = m7 -^ 
00 (in lattice units). For finite ^, however, we 
have artifacts of 0{a^). The Z{2) symmetry 
shields the system from 0{a) artifacts. 

We don't go through the corresponding formu- 
lae for the small field approximation here but just 
quote the result that the scaling artifacts remain 
of 0(0?) H]. There is a slight improvement on 
the 0{a?) level, but this is not a great progress. 

Finally we consider the dispersion relation for 
the FPA: 
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For each li a pole contributes to the correlation 
function and hence 
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This is the exact continuum energy spectrum. It 
is remarkable that it is obtained although ki G 
] — TT, tt]. Here we also observe the restoration of 
continuous rotational symmetry resp. Lorentz in- 
variance (which can not be seen directly in a{k).) 
Moreover also the continuous translation invari- 
ance is restored (note that fci + 27r/i covers the 
whole real axis). 



Using the perfect operator for the chiral con- 
densate (^l]) and inserting the FPA's gap eq. 
(|l|), we obtain also for the FPA ^ 
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Now the dimensionless ratio becomes 
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for any ^, i.e. we have perfect scaling. 

The result that in the FPA cut-off effects are 
completely eliminated means that the perfect 
classical actions are situated on the renormalized 
trajectory. This can be understood from the fact 
that in the limit N ^>- oo the path integral reduces 
to the classical solution; all fluctuations around 
the minimum are suppressed with weight zero. 

9. Summary and outlook 

What we should know to determine the clas- 
sically perfect action for any N is the FPA 
^S*[x,X,^], which is reproduced under RGT. 
However, this could be obtained analytically only 
in the small <& approximation. 

If we let N ^ oo, it suffices to consider a con- 
stant $0 and S* can be calculated exactly (non- 
perturbatively) . Then the large TV limit does a 
second job for us and makes the line of perfect 
classical actions coincide with the renormalized 
trajectory, such that we indeed observe perfect 
scaling. 

For finite N, both trajectories will change 
(there occur higher order interactions etc.) and 
we have no reason to expect that they still coin- 
cide. It would be interesting to determine the size 
of the cut-off artifacts on the classical trajectory 
in this case. This could be done either numer- 
ically for finite N, similar to [|l|, or by a 1/A^ 
expansion. A basis for the latter is given by our 
results for the small field approximation applied 
to the non-zero modes of $(fc). 

We are confident that our result holds (qualita- 
tively) also for 0{N) and CP{N) models. Then 
for 0{N) the success of the work of Hasenfratz 

^Both equations are quite complicated, but they have 
exactly the same structure, which allows for a drastic 
simplification. 



and Niedermayer indicates that the artifacts are 
still very small at A^ = 3. 

10. Conclusions 

At A^ — > oo the procedure of Hasenfratz and 
Niedermayer leads to a perfect scaling in the 
Gross-Neveu model, as we showed analytically. 
In the scaling behavior, cut-off artifacts are com- 
pletely eliminated at any correlation length. That 
the FPA really describes continuum physics on a 
lattice (of finite spacing) can also be seen from 
the fact that the continuum energy spectrum, 
and with it continuous Poincare invariance, are 
restored. 

The small field approximation does not im- 
prove the scaling essentially. Here the desired 
direction out of the critical surface is projected 
on a subspace, which is too restrictive. 

Also the asymptotic scaling is improved for the 
FPA, but this is not of primary importance. 
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